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ABSTRACT 



Tills paper investigates the feasibility of an air-to- 
surface missile delivery system employing impact prediction 
to assist a human operator in controlling the system. The 
impact prediction information was used to drive a display 
wnich the operator used in controlling the flight of the 
weapon. It was desired to find whether impact prediction 
proved useful in weapons control. 
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I . INTRODUCTION AND STATEMENT OF THE PROBLEM 

A. TYPICAL ASM WEAPONS SYSTEMS 

Most older air-to-surface weapons share the characteris- 
tic that once the weapon is launched it is under the influence 
of only the environment, e.g., gravity, wind drift, drag, etc. 
Examples in this category of weapons are bullets, general- 
purpose bombs, and rockets. This makes it the responsibility 
of the weapons officer to insure that the fire-control solu- 
tion is accurate enough for the weapon to impact at the tar- 
get within the kill radius of the warhead. In the normal 
tactical situation this is a difficult task to perform prop- 
erly. As a result some of the newer weapons systems have 
been designed with a capability for guidance of the weapon 
to the target. 

These newer systems require that to provide guidance for 
the weapon, the weapons control officer must maintain visual 
contact with the weapon during its flight. Again, this solu- 
tion to the problem suffers from tactical restrictions. In 
a combat environment it is very difficult to maintain a flight 
path which allows the weapons control officer to maintain 
visual contact with the weapon during its flight. 

In order to overcome the shortcomings of these systems, 
even newer systems such as the Phoenix employ very sophisti- 
cated missiles which allow the launching aircraft to stand-off 
from the target and still deliver a weapon so as to destroy 
or seriously damage the target. One aspect of these systems 
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which is objectionable is the high cost of the missiles which 
are employed. As a result this study was made of a system 
which would employ a stand-off delivery capability while em- 
ploying a relatively unsophisticated, and therefore inexpen- 
sive missile to deliver the warhead. 

B. ANALYSIS OF STAND-OFF WEAPONS DELIVERY 

Stand-off weapons delivery is best defined as the ability 
to remain outside the destructive range of the enemy's weapons 
and deliver ordnance on the enemy installation. The advan- 
tages of remaining outside the destructive envelope of the 
enemy's weapons are obvious. The disadvantages of standing 
off from the target are numerous. The greater the stand-off 
distance, the more accurate the guidance and system signals 
must be as well as the knowledge of the target location. 

This is necessary to provide a reasonable probability of tar- 
get kill. 

In its most basic form the air-to-surf ace weapons delivery 
problem consists of two parts: the detection and location of 

the target by some means, and delivery of a weapon to the 
target. 

The solution which is proposed in this paper is based on 
three assumptions about what the author feels an air-to- 
surf ace weapons system should incorporate. These three assump 
tions are 1) that the target should be located by some method 
which does not endanger the launching aircraft; 2) that the 
primary purpose of air-to-surf ace missiles is to deliver ord- 
nance on the target, not to carry a large amount of expensive 
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electronic equipment to be blown to pieces along with the 
target; and 3) that to improve the overall accuracy of the 
system, the missile will be under some form of command guid- 
ance during its flight to the target. 

The solution which is proposed in this paper is one of 
the "shoot and loiter" school of weapons control. Those of 
this line of thought contend that it is more economical to 
shoot a relatively unintelligent missile at the target and 
maintain the sophisticated system parts in the launching air- 
craft. This requires that the launching aircraft remain in 
the target area until the weapon has hit the target which 
increases the probability of the launching aircraft being 
destroyed. 

The opposite point of view is that the aircraft should 
launch the missile and get out of the target area. This 
approach necessitates the use of a very sophisticated missile 
but reduces the probability of the launching vehicle being 
destroyed. 

The trade-off to be considered between the two approaches 
is to obtain the proper balance of economy versus probability 
of survival of the launching aircraft. 
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II . THE PROBLEM; ANALYSIS AND SOLUTION 

In this chapter the two divisions of the problem will be 
discussed in detail. The first is the problem of target loca- 
tion. The second, delivery and guidance of the weapon, will 
be subdivided into two sections. The first section will be 
concerned with the initial fire-control solution for the mis- 
sile flight path and the second with the guidance of the 
missile to the target. 

A. LOCATING THE TARGET 

The task of locating the target will be limited to the 
case of a radiating target such as a radar. This limitation 
was made in order to propose a solution in which the launch- 
ing aircraft both locates and destroys the target. The 
assumption of a radiating target allows the use of a direction- 
finding method to locate the target. 

Several methods of direction finding are available. These 
include triangulation of radiation given off by the target 
using two or more receivers at known locations, the use of 
filtered sequential measurements of bearing and depression 
angles from a single aircraft at known positions, or even the 
simple method of having two ground observers locate the tar- 
get. For this solution the second method was used. 

1 . Basic Geometry of Target Location 

In Fig. 1 the basic principle of location by triangu- 
lation is illustrated. The points 1 and 2 could be either 
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Figure 1. Target Location by Triangulation 
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two distinct receivers listening to a signal at the same time, 
or they could be the two locations of a single aircraft at 
different times. 

Figure 2 gives a definition of the angles involved in 
all the target-location schemes. The azimuth angle is that 
angle measured from a true north reference positive clockwise 
to the line-of-sight (LOS) between the aircraft and the target. 
The depression angle is the angle measured from a reference 
plane normal to a line connecting the aircraft and the center 
of the earth. An angle which rotates below the normal plane 
is considered negative and one which rotates above the plane 
is considered positive. Appendix E contains a discussion of 
the reference systems used in the entire problem. 

In the following two sections two different methods 
of location of the target by use of angular information will 
be presented. 

2 . Location Using Positional and Angular Infomation 

The first of the two methods for target location will 
be explored here for the two-dimensional case to show the 
principle involved. An analysis of the three-dimensional 
case is contained in Appendix A. 

Figure 3 shows the geometry of the two-position method 
of target location. The aircraft, at two different times and 
at different positions, measures the azimuth angle from the 
aircraft to the target. At the same time as the angle is 
measured, the position of the aircraft is recorded. Assume 
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Figure 3. Two-Position, Planar Target Location 



The measurements 



that the target coordinates are and 
are listed below. 

Time Aircraft Location Azimuth Angle 



1 

2 




0 

0 



1 

2 



From Figure 3 it can be seen that 

Xf-Xi 

Tan 0, = - 

^ ^t ^1 



V^2 

Tan 0^ = - j - - 
^t 2 



(II-l) 



(II-2) 



Writing the two equations as simultaneous equations 
and solving them: 



\ ®1 
Xt - X 2 = (Y^-Y^) Tan 0^ 



X^ - Y^ Tan Tan 0^_ 



X^ - Y^ Tan ^2 ^ ^2 " ^2 ®2 * 



(II-3) 



(II-4) 



Letting A = Tan 0^^, B = Tan 02/ C = X^ - Y^^ Tan 0^^, 
and D = X 2 - Y 2 Tan 02 and writing the equations in matrix 



form 



1 

1 




(II-5) 



Using normal matrix algebra techniques, the solution 



can be shown to be 




DA-BC 

A-B 




D-C 

A-B 



(II-6) 

(II-7) 
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Expanding the results into their full form by replacing 



the A, B, C, and D with their respective values gives 



X 



t 



Tan0j^ -TanG^ (Xj^-Y^^TanG^) 

TanG^-TanG2 



(II-8) 



and 




(X2"Y2TanG2) - (X^^-Y^TanG^ ) 
TanG j^-TanG 2 



(II-9) 



This technique of locating the target can be quite 
effective and accurate if the two positions used in deriving 
the solution are far enough apart and if the bearing angle 
measurements are respectably accurate. Ideally the two bear- 
ing angles should be separated by ninety degrees. This sep- 
aration will reduce to a minimum the effects of any noise in 
the measurement of the angles. If the target radiates for 
only short periods of time and the ambient conditions are 
such that the measurements are very noisy, the derived target 
location could be greatly in error. For this reason an al- 
ternate approach to the problem of locating the target was 
taken. The alternate method uses only angular information 
and a hybrid filtering technique. 

3 . Location Using Angular Information and Hybrid Filtering 
The target location method discussed in this section 
is the result of research done by Demetry [1]. The method em- 
ploys an unusual hybrid filtering technique in which a Kalman 
filter is considered linear for calculation of a gain sched- 
ule but in which prediction is done nonlinearly. The hybrid 
filter will be discussed here briefly. 
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It is desired to locate a radiating target by a method 
in which only a bearing angle to the target is measurable. 

The equations below give the relationships used in the formu- 
lation of the filter. 



-k|k -^®k|k'®k|k'ik^ 



(II-IO) 



where 



a = 



XT 



k k 



YT 



k 1 k 



(II-ll) 



is the estimated target location at time k and 




■x(t)' 

X (t) 

y (t) 
y(t) 



(11-12) 



is the aircraft vector. 

The aircraft position and velocity components are 
considered to be measurable without noise for this evalua- 
tion as they are throughout the remainder of the problem. 

The bearing angle dynamics are described by the 
differential equations 



ip (t) 



0 ( t) 

5(t) 

6(t) 



f (a,C(t) ,t) 



(11-13) 



where ^ is a nonlinear vector function. The angle measurement 
at any discrete time k is 

^k " ^^k ’^k' ^ ^ ^ (11-14) 

in which H = [10 0] and Vj^ is the error in the bearing angle. 
The Vj^ is random in nature with an assumed mean value of zero. 

The gain schedule was calculated by assuming linear 
dynamics and constant acceleration for the process and using 
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the recursive Kalman gain equations [2] . However, the predic 
tion of angular information was done by using the actual anqu 
lar dynamics in a discretized form rather than employing the 
linear matrix which was used in computing the filter gains. 
These equations , by using true system dynamics to make approx 
imate predictions, reduce the location error since, in effect 
they consider the system linear. In this system this is not 
a bad assumption since the sample period is relatively short 
compared to the rates of change of the system variables. 



®k+l|k 



, I XT. I . -X, 
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(11-15) 
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"yk"‘%ik-^ki‘%ik-yk 

"^k|k-yk-yk^'"^'®k|k 



-Xkj2)2 



(11-17) 



This is a brief discussion of this particular method 
of location of a target. For the reader interested in pursu- 
ing this further. Ref. 1 will be of interest. 



B. THE FIRE-CONTROL SOLUTION 

Before it is possible to launch a missile at a target it 
is necessary to develop a solution for a flight path which 
will take the missile from the launch point to the target. 
The accuracy with which the solution need be computed is en- 
tirely a function of the amount of guidance which will be 
available to influence the missile during its flight. Also, 
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before an initial solution can be computed, it is necessary 
to know what approximate flight path the missile will follow 
from the launch point to the target. 

In the solution which is being investigated here, the 
proposed flight path for the missile is a nominal parabolic 
path. The path the missile actually travels will not be para- 
bolic due to air drag, control effects, and other outside 
influences . 

There are two main reasons for selecting a parabola as a 
nominal flight path. The first concerns a relationship be- 
tween the height of a parabola and distance to the target 
from a point under the peak of the parabola. This ratio is 
referred to as the A/R ratio and its importance is discussed 
in Appendix B. The second consideration is that the missile 
will be under the influence of gravity in a very helpful man- 
ner. The gravitational force will help to maintain the veloc- 
ity of the missile and hence, the controllability of the 
missile, after the boosted phase of the missile's flight. 

1 . Geometry of the Fire Control Solution 

The general geometry of the fire control problem is 

illustrated in Fig. 4. In the solution the aircraft, point 

A/C on the figure, is considered to be at the X origin of the 

figure and the target, point T in the figure, is considered to 

lie on the X-axis. The parabola will be defined such that it 

contains the two points of the launching aircraft (A/C) and 

the target (T) . Since a parabola is defined by a second-order 

2 

polynomial, e.g., of the form y = ax + bx + c, some third 
point on the parabola must be defined. 
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A 

In Fig. 5, the parabola of Fig. 4 is redrawn, this 
time witn three points and four distances defined. For the 
purpose of defining the nominal parabola, the launching air- 
craft is placed on the Z axis so that range is zero at the 
origin. The altitude of the aircraft is labeled zA. A sec- 
ond point is the location of the target. The target is at 
zero elevation and a distance D from the origin. The third 
point, which is necessary for the parabola to be defined, is 
the peak of the parabola. The altitude of this point is A 
and it is D - R from the origin. Using these three points 
it is possible to solve for the coefficients of the second- 
order polynomial defining a parabola. 

Writing the equations of the parabola at the three 
points gives 



a(0) ^ 


+ b(0) + 


c = 


zA 


(11-18) 


a(D-R) ^ 


+ b(D-R) 


+ c 


= A 


(11-19) 


a(D) ^ 


+ b(D) + 


c = 


0 


(11-20) 




R = kD 


• 




(11-21) 



Equation 11-21 is used to define the relationship 
between the range from the origin to the target and the 
range from under the peak of the parabola to the target. 

2 . Method of Obtaining the Launch Angle 

In a system in which the launch point and target point 
are coplanar in the earth plane, a ballistic mass will travel 
furthest when the velocity vector is initially directed 45 
degrees above the horizontal. This. is true as long as the 
mass has no external forces acting on it other than gravity. 



20 



Two factors involved in this investigation make the 
choice of a 45-degree launch angle unacceptable. The first 
is tiiat the launching point and target are not earth coplanar 
and, secondly, that the mass is subject to thrust from a 
rocket motor. Thus, the problem of computing the launch angle 
becomes one of maximizing the range as a function of the dis- 
tance between the target and launching point in the Z plane. 

The results of this particular maximization were 
interesting. For aircraft locations which were up to 5000 
meters noncoplanar with the target, the launch angle for max- 
imum distance was found to vary only one degree, between 47 
and 48 degrees. (The launch angle as used here is defined as 
the angle between the velocity vector and the horizontal at 
the launch point. This is not the angle of the velocity vec- 
tor at the time of burnout which is the angle which deter- 
mines the maximum distance of flight.) The launch angle was 
set at a constant 47.5 degrees as a result. This angle should 
give the maximum range regardless of range to the target. 
Starting with the launch angle that gives maximum range should 
prove entirely acceptable since one of the assumptions about 
the system was that the missile would be under some type of 
guidance during its flight. 

C. GUIDANCE BY THE HUMAN OPERATOR 

Tnroughout this paper mention has been made of the fact 
that the missile will be guided by the human operator through- 
out its flight. In actuality the missile was not released to 
the operator's control until after it had finished the boosted 
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phase of its flight. During this time the system turned the 
missile up to the launch angle and yawed it to the proper 
flight heading. Once the missile started its free flight, 
a display of the predicted impact point of the missile rela- 
tive to the target was generated for the missile control 
officer (MCO) . At the same time, the control surfaces of the 
missile were set to a neutral position and control of the 
missile was given to the MCO who then guided the missile the 
remainder of the distance to the target. 

In order for the predicted missile impact point to be 
displayed, it was first necessary to track the missile and 
based on the tracking information, predict its impact point. 
The prediction information was then presented to the MCO in 
such a manner that he could make decisions as to the guidance 
control to be given to the missile. 

1 . Tracking the Missile in Flight 

In this problem the missile was tracked in flight by 
a radar which was capable of providing four items of informa- 
tion. These were the range to the missile, the range rate of 
change to the missile, the azimuth angle to the missile, and 
tne depression angle to the missile. The reference system by 
whicn the angles were measured is discussed in Appendix E. 

The range was measured in the positive sense while the range 
rate was considered positive for a missile whose range was 
increasing and negative for a missile which was closing the 
launching aircraft. 
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since the radar was operating in a natural environ- 
ment, the information provided by the radar was not perfect. 
This noise would introduce errors in the estimation of the 
location of the missile during its flight, so it was desirable 
to use some technique to reduce the effect of the noise on 
the measurements. The technique used was Kalman filtering. 
Before deriving the necessary Kalman filter for this problem, 
the desired missile information vector should be investigated. 

Figure 6 presents the geometry for developing the 
missile information vector from the spherical vector of the 
radar. The desired missile information vector contains six 
components: the XYZ cartesian location of the missile rela- 

tive to the aircraft and the XYZ cartesian velocities, again 
relative to the aircraft. The acceptability of providing 
only a six-element state vector for the missile is dependent 
upon the sample rate at which the state vector is developed. 
Even though the missile is being acted upon by external forces, 
it is possible to consider the missile to be a piece-wise 
constant-velocity missile in either spherical or cartesian 
coordinates if the interval between samples is sufficiently 
small . 

From Fig. 6 the equations for the missile location 
relative to the aircraft can be developed. Then by the use 
of the time derivatives of the equations of location, the 
equations of the cartesian velocities can be found. First, 
the equations of location. 
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Figure 6. Geometry of Missile Location 



R sin (}) 



( 11 - 22 ) 



Z = 
m 

X = R' sin 6 = R cos <P sin 6 (11-23) 

m 

Y = R' cos e = R cos (p cos e (11-24) 

m ^ 

For the case shown in Fig. 6, the angle <p is negative, 
giving the proper negative value for . Likewise, the signs 
of the remaining two values will be correct depending on the 
quadrant of the angle 9. 

To form the equations of the velocities in the carte- 
sian frame , it is necessary only to take the time derivatives 

of equations 11-22, 11-23, and 11-24. These are given below. 

• • • 

Z = R COS d) (|) - R sin d> (11-25) 

m 

• • • • 

X = R cos d) cos 9 6 - R sin d) sin 9 d) + R cos d) sin 9 
m 

( 11 - 26 ) 

• • • • 

Y = -R COS d) sin 9 9 - R sin d) cos 9 d) + R cos d) cos 9 

r -t- r r 

(11-27) 

From equations 11-22 through 11-27 a list of the 
necessary spherical values of missile movement to develop 
the cartesian state vector can be collected. These equations 
show that it is necessary to have a vector of spherical values 
of range, azimuth, depression, and the associated rates of 
change of the three quantities. From this requirement a Kal- 
man filter can be developed which will not only reduce the 
effect of the measurement noise on the quantities which are 
being measured, but also provide estimates of the angular 
rates which are not measured. 

The Kalman filter for filtering and estimating the 
radar vector has the general Kalman predictor-corrector form 
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of equation 11-28. The Gj^ is computed from the recursive 
equations of 11-29, 11-30, and 11-31. More complete details 
of the derivation of the filter and the values used in the 
gain schedule computation are given in Appendix C. 





(11-28) 




(11-29) 


k|k " ^^“\^^^k|k-l 


(11-30) 


k+l|k " ^^klk*^ ^ 


(11-31) 



2 . Predicting the Missile Impact Point 

In predicting the missile impact point it is desirable 
to have an algorithm which is fairly accurate. This is the 
case since it is the point at which the missile is predicted 
to strike that will be used by the MCO to develop his guidance 
commands. A requisite increase in accuracy over a ballistic 
model can be obtained by introducing the effects of drag into 
the equations of motion for the missile. For the prediction 
algorithm developed here, the drag was considered to be oro- 
portional to the velocity of the missile. 

In this solution the drag coefficient of the missile 
is taken to be constant from the time of the prediction until 
the missile strikes the target. This, of course, is an invalid 
assumption but one which is necessary in order to develop a 
closed- form solution for the prediction equations. The error 
which is introduced by this assumption will have the effect 
of causing the MCO to hold the missile aloft longer than would 
normally be necessary. This action by the MCO will tend to 
increase the A/R ratio, thus making the missile impact on the 



26 



target more vertically and hence less subject to error from 
target altitude location error. (See Appendix B) . 

a. Equations of Motion for Impact Prediction 

The equations of motion for impact prediction were 
written for each of the three component directions of the 
cartesian reference system. The equations were based on an 
assumed point-mass model of the missile subjected only to 
drag and gravity forces. The X and Y equations are similar 
to each other, varying only in the drag coefficient. The Z 
equation takes into account the effect of gravity. The equa- 
tions of motion are 



mX = -K^X 



(11-32) 



where 



my = -K 2 Y 



mZ = -K^Z - mG 



= p|x|a^ 

= p I Y |A 
K3 - p|z|a^ 



(11-33) 

(11-34) 

(11-35) 

(11-36) 

(11-37) 



and m = missile mass, G ^ gravity, and A^ = reference area. 

The method for predicting the impact point of the 
missile is to first compute the drag coefficients (K^^, , 

and K^) which are then considered constant until impact. The 
Z motion equation is then solved for the time when the missile 
Z location is equal to the target Z. First, it is necessary 
to solve equation 11-34 for the missile Z as a function of 
time . 
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b. 



Impact Prediction Equations 

The result of solving for the Z motion as a func- 



tion of time is 



Z (t) = (1-A) e - Bt + C 



where 



A = 



mZ (0) ^ m^G 



K. 



B = 



mG 

iTT 



c = Z (0) 
m 



a = 



K. 



m 



(11-38) 

(11-39) 

(11-40) 

(11-41) 

(11-42) 



This mixed equation is not directly solvable for 

that value of time such that Z is equal to Z . . The Newton- 

m ^ t 

Raphson iterative method is employed to find the value of 

time which will satisfy the equation for a specific set of 

conditions (a particular Z (0) and Z (0) ) . The value of time 

mm 

thus found is substituted into the X and Y solutions to find 
the predicted impact point of the missile. The X and Y 
predicted positions at impact are given by equations 11-43 
and 11-44. 



mX (0) 



X = X (0) + ^ 

p m K 



m 



-Kt t. /m 

/T 1 1 \ 

(1-e 



mY (0) 



Y = Y (0) + „ 

p m K 



m 



(1-e 



-K^ t^/m 



(11-43) 



(11-44) 



where t^ is the value of impact time from the solution of 
equation 11-38. 
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3. 



Display of the Situation 



In order for the MCO to make intelligent decisions 
regarding the application of control to the missile during 
flight, there had to be some method for displaying the pre- 
dicted results of the current flight path. What was displayed 
to the MCO was the current prediction of the impact point 
of the missile assuming that there would be no further control 
exercised over the missile. The display is generated such 
that the center of the screen corresponds to the last esti- 
mate of the target position and the predicted impact point is 
displayed relative to the center of the screen. The disolay 
is also rotated so that motion of the impact point will 
correspond to control stick motion. Figure 7 illustrates 
the display screen and Appendix H further discusses the man/ 
system interface. 
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Figure 7. The Situation Display 



Ill . 



THE SIMULATION OF THE PROBLEM 



A proposed solution to a problem cannot be considered to 
be valid unless it is tested in some way in practice. One 
method of evaluation is to build an actual system and test 
it in a practical operating environment. This is generally 
too expensive to be feasible. Another method of testing a 
proposal is to simulate the system on an appropriate computer. 
The second method was chosen to evaluate the problem solution 
proposed in this paper. 

This simulation was done on two different computers. The 
first was an IBM System 360 Model 67, operating under the MVT 
operating system. This system was used for initial program 
debugging since the Operating System has an excellent error- 
monitoring system. The second system was a hybrid computer 
consisting of an SDS 9300 digital computer and a COMCOR CI- 
5000 analog computer. 

A. SIMULATION OF TARGET LOCATION 

The simulation of the target location was accomplished 
by using the results of work in passive target location done 
by Demetry [1]. The results of the hybrid location method 
discussed in Chapter II were used by making use of the calcu- 
lated mean-square error values in target X and Y location 
as a function of time, from the hybrid location algorithm. 

Using these values involved a process of both normaliza- 
tion and interpolation. Since the simulation is designed for 
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a time increment of one second, the target location mean- 
square errors had to be interpolated since they were computed 
on a time increment of two seconds between , observations . 

Tne interpolation was made on the assumption that the 
mean-square errors of location were piece-wise linear over 
tne measurement interval involved. The normalization process 
was necessary since when the values were generated, the air- 
craft was flying a set course which had a direct effect on 
tne value of the mean-square error. Thus, the values were 
normalized by multiplying them by a factor which was the 
closest point of approach (CPA) to the target for the prob- 
lem divided by the CPA to the target in the mean-square-error 
computation program. 

B. SIMULATION OF THE FIRE-CONTROL SOLUTION 

A real-time simulation of the fire-control solution was 
eliminated after some research into the problem. A constant 
launch angle was used in its place. The solution was arrived 
at Dy assuming the missile to be a drag-free ballistic mass 
with the propulsion system of the true missile attached. 

Tnis approximation was used in an iterative method of comput- 
ing the angle which would give the maximum flight down range 
for a particular launch altitude. The iteration was performed 
on both the launch altitude and the launch angle. For this 
investigation the angular increment was taken as one degree 
and the altitude increment was fifty meters. 

For each launch altitude, the angle which gave the maximum 
down-range flight of the assumed missile was plotted against 
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the launch altitude. The results were surprising in that it 
was found that for altitudes between fifty and five thousand 
meters, the angle varied only between 47 and 48 degrees. The 
launcn angle is defined as the initial angle between the hor- 
izontal and the velocity vector. This initial angle is re- 
duced by the effect of gravity to an angle at burnout which 
will give the maximum down-range flight of the mass. The 
force of gravity is thus responsible for the initial angle 
being greater than 45 degrees. 

It must be noted that this angle was chosen to give the 
maximum down-range flight of the assumed missile. Thus, by 
shooting the true missile at this angle the maximum possible 
down-range flight would be possible for zero control. By 
the application of control to the missile in the pitch mode, 
this range can be either increased or decreased. In the 
desired situation this range should have to be decreased by 
using down pitch. This action would cause an increase in 
tne A/R ratio which has already been demonstrated to be 
desirable . 

C. SIMULATION OF THE MISSILE 

The missile model was the most important part of the simu- 
lation and thus was treated with the most rigor. The attempt 
was made to use a fairly accurate model in that dynamics 
approximating an actual missile were assumed. Appendix F 
contains a description of the missile model. 
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D. GUIDANCE SYNTHESIS 



The guidance synthesis for the simulation was done in two 
parts : one part for the simulation while it was executed on 

i" 

tne System 360 and the second for the actual data-gathering 
simulation on the SDS/COMCOR hybrid. 

1 . System 360 Guidance Synthesis 

The guidance synthesis used during the preliminary 
runs on the System 360 was a primitive, bang-bang control. 

Yaw was either plus or minus two degrees per second and pitch 
was a plus or minus 0.5 radian deflection of the control sur- 
face. The assumptions on which this control was synthesized 
were the same as those discussed in the following sections. 

2 . SDS/COMCOR Guidance Synthesis 

Before going into the discussion of the methods used 
to develop the guidance commands in the hybrid simulation, the 
assumptions on which the command synthesis was based will be 
discussed. 

a. Roll-Independent Control 

The assumption that control is roll independent 
is valid since the majority of missile control systems operate 
in this manner. Also the assumption allows a simplification 
of the command vector synthesis. The missile is thus consid- 
ered to either be roll stabilized or the control system in 
the missile is so designed to compensate for any roll of the 
missile . 
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b. Two-Component Control Vector 

With roll independence assumed, the control vector 
consists of two components, a pitch and a yaw component. 

These components were limited in their magnitude by the syn- 
thesis program. 

The yaw was limited to a maximum rate of twenty 
degrees per second by the control mechanism itself. The 
pitch was limited to plus or minus forty-five degrees of 
control deflection. Additional limits on the control deflec- 
tions were built into the missile model. Yaw rate was lim- 
ited so that radial acceleration on the missile would not 
exceed ten G's. 

c. Control Vector Formation 

The control vector was formed by the human opera- 
tor directing a control stick. The maximum deflection of the 
stick corresponds to the maximum amount of control which 
could be generated. This control was then applied to the 
missile each second until the control stick was neutralized. 
The control affected the missile in that the control surfaces 
of the missile remained deflected only as long as the stick 
was held in a non-neutral position. 
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IV. DISCUSSION 



The simulation of this problem was performed only a few 
times with the completed program. As a result it was not 
possible to draw concrete conclusions about the validity of 
the total concept tested. However, several observations were 
made about the program and some areas where further investi- 
gation could be made were discovered. 

A. OBSERVATIONS 

One observation from the running of the program supported 
tne supposition that the use of impact prediction in an air- 
to-surface environment was effective. The errors in impact 
prediction seemed to be related closely to errors in producing 
the measured missile state vector. Overall the impact pre- 
diction approach to missile control appeared very feasible. 

Another observation which was made concerned the initial 
laxinch angle assumed for the simulation. Apparently the use 
of a fixed angle of launch intended to give maximum down- 
range flight was not as good an . assumption as first thought. 
For targets which were close to the launching aircraft, it 
was possible to overfly the target by shooting for maximum 
range. It would be much more desirable to have the launch 
angle be a function of the initial range to the target. 

The overfly or fly-by effect can be a serious problem. 

If the missile bypasses the target with a small angular offset 
to either side, it is possible that immediate application of 
yaw control in such circumstances will cause the missile to 
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circle the target at a constant range. This can occur since 
the missile is limited in the amount of radial acceleration 
it can withstand. It should be possible to overcome the prob 
lem by the use of tactical doctrine if the missile does over- 
fly the target. One suggested doctrine would be to continue 
the fly-by course until the missile is far enough out to make 
the turn and hit the target area. This maneuver is illus- 
trated in Fig. 8. 

B. FURTHER INVESTIGATION 

Three areas exist which the author feels merit further 
investigation. The areas are 1) prediction accuracy and 
methods; 2) radar filter accuracy; and 3) display methods. 

The accuracy of impact prediction should be investigated 
further. Although it was shown that the impact prediction 
method worked, no study of the accuracy of the method em- 
ployed was done. Also, it is possible that different meth- 
ods of impact prediction would prove to be more accurate 
than the velocity-proportional model used in the simulation. 

Prediction accuracy, from all indications, was very 
dependent on the accuracy with which the missile state was 
estimated. The key portion of the program in the missile- 
state estimation was the radar vector filter. The filter 
used in the problem was based on the assumption of a constant 
velocity model of the missile. The use of a filter based on 
constant acceleration should be investigated as lags were 
noted in the measured missile state compared to the true mis- 
sile vector when the missile was making hard or continued 



maneuvers . 
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The method of displaying the prediction information to 
the MCO also deserves further investigation. The display as 
it was implemented did not provide any indication of the veloc- 
ity of the missile. It is felt that if the MCO were given 
missile velocity information, he would be able to better con- 
trol the flight of the missile. This is especially true in 
the fly-by case. 

As the missile bypasses the target, the display gives a 
false indication as to the control necessary to hit the tar- 
get. This problem is generated as a result of the almost 
180-degree change in angle by which the display is rotated. 

Figure 9 illustrates an improved display which should 
solve the rotation problems caused by fly-by and also give 
velocity vector information to the MCO. The display is of 
the entire area of the problem. The axes of the new display 
are neither translated nor rotated as they were in the 
simulation. 

Three points are included on the display. One is the 
current XY measured position of the missile. This is the 
left-most point in Fig. 9. The other two points are the 
current predicted impact point and the last estimate of the 
target location. A line from the missile position to the im- 
pact point will indicate the velocity information desired. 

Also the non-rotation of the display axes should alleviate 
the false indication of necessary command. 
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Figure 8. Course for Fly-by Error Correction 




Figure 9. Suggested Improved Display 
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C. RESULTS OF A RUN 

The figures in this section show the results of only one 
run of the simulation program. The figures give two planar 
plots of the flight path. One is a plot of the XY plane 
flight of the missile (Fig. 10) and the second is a vertical 
plane view of the flight path (Fig. 11) . 

Other plots are included which show the control vector 
as a function of time (Fig. 12) and the missile location 
errors as a function of time (Fig. 13-14) . 

It should be emphasized again that this is only one run 
and should not be considered to be in any way conclusive. To 
reach any definite conclusions a much more exhaustive investi- 
gation would have to be made. 
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Figure 10. Representative Run XY Flight Profile 
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Figure 11. Representative Run RZ Flight Profile 
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Figure 12. Missile Control vs. Flight Time 
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Figure 13. Missile Position Error vs. Flight Time 
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Figure 14. Missile Velocity Error vs. Flight Time 



APPENDIX A 



TARGET LOCATION BY SUCCESSIVE ANGULAR MEASUREMENTS 
FOR THE THREE-DIMENSIONAL CASE 



In chapter II there was given a brief description of 
locating a target in the two-dimensional plane by processing 
two sets of angular measurements and two associated aircraft 
positions. In this appendix the two-dimensional case will 
be expanded to the third dimension. 

Figure 15 shows the three-dimensional geometry involved 
in this problem. This figure is then shown in two plane views 
the XY plane projection in Fig. 16 and the RZ plane projec- 
tion of Fig. 17. 

From Fig. 16 the two locations of the launching aircraft 
are marked as X^^ , and X 2 , '^2 '^ith the associated azimuth 
angles 0^^ and 62 . Using the geometric relationship for tan- 
gent, the tangents of the angles shown are 



v^_y^ = tan 0 
t 1 




tan 02 



(A-1) 



(A-2) 



Treating the equations as simultaneous equations and setting 
them up for solution, 



- Y^ tan = Xj^ - Y^^ tan 0^ (A-3) 

X^ - Y^ tan 02 ^ ^2 " ^2 ®2 * (A-4) 
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Letting tan 6^^ = A, tan &2 = B, tan 6^ = C, and 

X 2 - Y 2 tan © 2=0 and substituting these values into equa- 
tions A-3 and A-4 results in 



X, - AY. = C 
t t 



X^ - by. = D . 
t t 



(A-5) 

(A-6) 



Rewriting these equations in matrix form gives 



'1 


-A' 
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'c‘ 
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-B 


!y^ 
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- t- 







(A-7) 



which has a solution of the form 
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-A ' 'C' 






! = 




i 
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(A-8) 



The solution of this equation, arrived at by normal matrix 
algebra, is 



(A-9) 
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i-B 
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’C‘ 




A-B 
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-1 
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L J 



Making substitutions back into the solution, the X^ and Y^ 
of the target are 

tan 6 (X„-Y„tan9«) - tan 0„ (X,-Y,tan6,) 

x^ = i r (A-10) 

t tan 6, - tan 






(X2-Y2tan82) - (x.-y tane.) 

tan - tan ©2 



(A-11) 



The solution for locating the Z coordinate of the target 
will now be derived. 

Figure 17 will be used to illustrate the derivation of 
the target Z location problem. Again using the tangent 
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Figure 15. Three-Dimensional Location Geometry 
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Figure 16. XY Plane View Location Geometry 
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relationships for the angles involved, it can be seen that 




tan 



(A-12) 




tan <t >2 



(A-13) 



From Fig. 16 the values for and can be seen to be 



The fact is very evident from equations A-14 and A-15 
that once the X and Y coordinates of the target have been 
found, any single measurement of the depression angle may 
be used along with the Z location of the aircraft at the 
time of the depression angle measurement to determine the Z 
coordinate of the target. 






(A-14) 



“2 = 



(A-15) 
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APPENDIX B 



EFFECT OF A/R RATIO AND TARGET Z LOCATION ERROR 
ON RANGE MISS DISTANCE 

Several references have been made in this paper to the 
A/R ratio and its importance in determining missile miss dis- 
tances at the target. This appendix is devoted to a discus- 
sion of the importance of this effect. 

A. DEFINITION OF A/R RATIO AND ITS EFFECT ON MISS DISTANCE 

The A/R ratio, as defined for a parabola, is the ratio of 
the height of the parabola at its apex to the distance from 
a point under the apex to one of the terminal points of the 
parabola. This ratio is illustrated in Fig. 18. It can be 
seen from the figure that the larger the A/R ratio, the more 
nearly perpendicular the terminus of the parabola at its 
intersection with the axis. The more perpendicular the ter- 
minus of the parabola, the less the point of intersection 
with the target plane will be affected by an error in locating 
the target plane. Figures 19 and 20 illustrate the relation- 
ship just described. 

In Fig. 19, which illustrates a small A/R ratio, it can 
be seen that for an error in target Z location, AZ , a rela- 
tively large error in range, AR, is generated. The alternate 
and desirable case of the large A/R ratio is shown in Fig. 20. 
For the same error in target Z location as in Fig. 19, the 
range error, AR, is much smaller than in the small A/R ratio 
example . 
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Figure 18. Illustration of A/R Ratio 
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Figure 19. Small A/R Ratio 




Figure 20. Large A/R Ratio 
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B. CURVES RELATING MISS DISTANCE AND A/R RATIO AS A FUNCTION 
OF TARGET AZ. 

The computation of a family of curves relating miss dis- 
tances and A/R ratio as a function of target AZ is straight- 
forward. Assuming a ballistic, drag-free mass, the time for 
that mass to fall the distance A under the influence of 
gravity was determined. This was done for the initial condi- 
tion of the mass having a zero Z velocity. In equation form 




2 

where Z = -9.8 m/sec and A is the initial Z of the mass with 
the target Z assumed to be zero. By using this time and the 
knowledge of the distance R from the apex of the parabola to 
the target, the initial R velocity was computed as 

. (B-2) 

1 

This velocity was then set as constant throughout the 
current run in which the Z location of the target was varied 
away from the perfect zero location for which the R velocity 
was computed. For each value of AZ , a time to fall, t^, was 
computed and then substituted into the equation 

R = V t. . (B-3) 

r f 

The difference between this distance thus computed and the R 
to the target was the miss distance. The series of computa- 
tions was repeated for a number of different A/R ratios and 
AZ 's. 

The results were then plotted with A/R ratio on the X 
scale and miss distance on the Y scale. The family of curves 
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of Fig. 21 was generated by connecting the miss distances at 
the different A/R ratios which were the result of the same 
AZ . 

The values used in this problem were 
R = 5000 meters (constant) 

A's such that A/R ratios of .1, .2, .3, .4, 

.5, .6, .7, .8, . 9 , 1.0, 1.5, 2.0, 2.5, 3.0, 

3.5, and 4.0 were generated. 

AZ's used were 1, 2, 5, 10, 20, 50, 100, 

200, and 500 meters. 
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Figure 21. Range Miss .Distance as a Function of 
AZ and A/R Ratio 



56 



APPENDIX C 



DERIVATION OF THE KALMAN FILTER 
FOR THE MISSILE -TRACKING RADAR [2] 

In chapter II the Kalman filter for the system radar was 
discussed. This appendix gives the assumptions which were 
used in the derivation of the filter. 

The filter was designed with the assumption that the 
missile-tracking radar provides measurements of range, range 
rate, azimuth angle, and depression angle to the missile 
relative to the aircraft. Another assumption was that the 
measurements were contaminated by additive noise whose sta- 
tistics were known. It was also assumed that the missile was 
not subjected to any random perturbations. Based on these 
assumptions, some of the matrices of equations C-1, C-3, and 
C-5 may be defined. 







(C-1) 


^k| k 


= '^-°k«>''klk-l 


(C-2) 


^k+l|k 


= H- Q 


(c-3) 


|X> 


= 4|k-i ^ °k<5.k-«ik|k-i> 


(C-4) 


1 — \ 
\ 

<x 


1 — 1 

1 

1 — 1 
1 

<xl 

o 

II 


(C-5) 



The Q matrices are identically zero since it was assumed 
that the missile was subject to no random perturbations. The 

^The symbols used in the equations C-1 through C-5 and 
elsewhere in this appendix differ somewhat from the svmbol- 
ogy used in R. C. K. Lee. [Ref. 2]. 
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R matrices are not zero since they convey to the filter infor- 
mation regarding the variance of the additive noise in the 
measured values. 

The 4> matrices are the same for all the filters since the 
same dynamics (constant velocity) were assumed for all the 
filters. Equation C-6 gives these matrices as 

n T 



$ = 4 ) = T 

rg az dp 



(C-6) 



The H matrices in equations C-1 , C-2 , and C-4 are measure- 
ment matrices which determine which variables are measurable 
by the filters. In the range filter is a 2 X 2 identity 

matrix. The H and H, are 1X2 matrices with the (1.1) 

az dp 

elements being unity. 

It should be noted that equations C-1, C-2, and C-3 are 
data-independent as long as the matrices involved in computing 
them are data-independent and if the H, R, and Q matrices 
are constants. If this is true, then the gain schedule can 
be computed offline from the actual simulation of the problem 
or operational environment and simply inserted into the com- 
puter memory for use at the proper time. With the precomputed 
gain schedule in use, the implementation of the filter to 
process data in real-time is a relatively easy matter and 
does not contribute greatly to the computational load on the 
computer. 

An important consideration in designing any filter is the 
expected value of the measured parameters which will be used 
to initialize the filter. This is especially important when 
choosing values for the initial covariance matrices, Pq_Iq- 
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As these filters were to be employed in conjunction with a 
missile-tracking radar, the assumptions were made that the 
initial range to the missile and range rate of the missile 
would be unknown and that the initial angles to the missile 
would also be unknown. These assumptions were reflected in 
the covariance matrix initializations when the gain schedule 
was computed. The values used are given for the range and 
range rate filter in equation C-7 and the angular filters in 
equation C-8. 



10 



10 



10000. m" 
0 

,2 



0 

1000. m^/sec^ 



0 



(C-7) 



10 rad^/sec^ 



(C-8) 



10 rad" 

0 

The values used for the R matrices in the gain schedule 
computation were 



R 



rg 



100 m' 
0 



0 

10 m^/sec^ 



R = R, = .0004 rad" 
az dp 



(C-9) 



(C-10) 



By using very large initial covariance matrices and set- 
ting the initial estimates of the state vectors to zero, the 
gain element (1,1) for the first time increment will be very 
close to unity, thus causing the first measurement of the 
particular state variable to become the first estimate of 
that variable. As time passes, however, the incoming meas- 
urement is weighted less and less as the gain goes down, i.e., 
the filter becomes more confident in the accuracy of its 
estimates . 
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However, this reduction of gain as a function of time can 
sometimes be self-defeating. In the case in which this filter 
was employed, the reduction of the gain with time could be 
disastrous. Some corrective action had to be taken to prevent 
the gains from approaching zero and, as a result, having the 
filter ignore new data. Otherwise, as the missile maneuvered, 
the filters would have very quickly begun to produce very bad 
estimates of the radar state vector. 

Thus, the effect which is desired is one which will not 
allow the gains to approach zero. There are two actions which 
will prevent this from occurring: letting the 0 matrices be 

non-zero, or truncating the gain schedule after some length 
of time kT . Since the Q matrices have already been considered 
to be zero, there remains gain-schedule truncation. By trun- 
cating the gain schedule the filters continue to accept in- 
coming measurements as containing significant information, 
and the amount of core storage in the computer is reduced 
greatly . 

The gain schedules in this particular simulation were 
truncated in such a manner that the previous estimate produced 
by the filters were weighted by a factor of approximately 0.8 
and the incoming measurements were weighted by a factor of 
approximately 0.2. 
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APPENDIX D 



DERIVATION OF THE SOLUTION TO THE 
IMPACT PREDICTION EQUATIONS 



This appendix discusses the derivation of the impact 
prediction equations used for the missile impact point pre- 
diction. Two methods are discussed: one in which the drag 
force is proportional to the velocity and a second in which 
the drag force is proportional to the velocity squared. The 
first was used in the simulation program. 

In order to derive closed-form solutions to the equa- 
tions of motion for the missile, it was assumed that the drag 
coefficient was constant over the prediction interval. 

This assumption will introduce errors into the predictions. 
However, this error should be reduced the closer the missile 
approaches to its terminal point. Near the terminal point of 
the flight the greatest error will be in the predicted time 
until impact since the Z acceleration will be the greatest 
at that time. This should not greatly affect the X and Y 
prediction since the X and Y accelerations will tend to be 
very small at that point. 

The method used to solve the equations is to first solve 
the Z motion equation for the time at which the missile Z is 
equal to the target Z. The value of time thus found, which 
is the predicted remaining flight time of the missile, is 
then used with the current X and Y positions and velocities 
to compute the predicted impact point. 
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A. METHOD I 



DRAG PROPORTIONAL TO VELOCITY 



In method I the drag force on the missile is considered 
to be proportional to the velocity of the missile. The equa- 
tions for this case are linear differential equations. 

The equations of motion for method I are 

• • • 

mX = -K^X (D-1) 

• • # 

mY = -K^Y (D-2) 

• • • 

mZ = -K^Z - mG (D-3) 



where 

K, = pIxIa (D-4) 

1 ' r 

* 

K2 = pIyIa^ (D-5) 

= p I Z I A (D-6) 

3 K I I r 

• • • 

in which the X, Y, and Z are the current values of the 

velocities of the missile, p is the atmospheric density, and 

2 

A is a reference area of .114 meter . 
r 

The solution for the Z motion equation D-3 has the form 



where 



0 = -Ae““^ - Bt + C 



A = 



mZ^(O) ^2^ 
m m G 



K. 



K. 



B = 



mG 



C = A + Z (0) - Z, 
m t 



a = 



K. 



m 



(D-7) 

(D-8) 

(D-9) 

(D-10) 

(D-11) 
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The problem now becomes the solution of the time-domain 
equation of the Z motion for the time at which the missile Z 
is equal to the target Z. The C coefficient in equation D-10 
is set to reflect this desired equality already. 

The method chosen to solve this equation is the Newton- 
Raphson method. Upon commencement of the iterations all of 
the required coefficients were computed using the current 
values of the positions and velocities. The termination test 
constant used to terminate the iterative process was 0.1 sec- 
ond. Whenever one value of time had changed from the previous 
value of time by an amount less than or equal to the test 
increment, the iterative process was terminated. The value 
of 0.1 seconds was chosen for the test value since regardless 
of the sample interval given to the simulation program the 
missile model was advanced in increments of 0.1 seconds. 

Thus, whenever the time was less than or equal to the minimum 
increment that the missile could be advanced, the iteration 
was stopped. 

The value of time which was returned from the Newton- 
Raphson iteration was then substituted into the equations for 
X and Y impact point prediction. These equations are 

X = X (0) + — ^ (1.0-e (D-12) 

pm K ^ 

« 

mY ( 0 ) 

Y = Y^(0) + — ^ d.O-e"^"^) (D-13) 

pm ix 2 

where 



3 




m 



(D-14) 



63 



Y 



(D-15) 




m 

These are the equations for method I prediction. The 

f' 

method II equations are more complicated to solve since they 
are non-linear. 



B. METHOD II - DRAG PROPORTIONAL TO VELOCITY SQUARED 

In this derivation of impact prediction equations the 
drag force on the missile is taken to be proportional to the 
square of the missile velocity. Again, the drag coefficients 
are computed at the beginning of the solution and considered 
constant throughout. 

The equations of motion are 

mic = -K^X^ (D-16) • 

• • * O 

mY = -K 2 Y (D-17) 

mZ = -K^Z^ - mG . (D-18) 

The steps of this solution are the same as used in the 
previous method: 1) solve the Z motion equation for Z as a 

function of time, 2) solve for the time when missile and 
target Z are equal, and 3) use the time thus found to compute 

X and Y from the prediction equations . 

P P 

First solve for Z as a function of time, writing equation 
D-18 as 



where 



and 



^ = 
dt 



-(Cy^+G) 




y 



= z 



(D-19) 



(D-20) 

(D-21) 
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Separating variables and integrating, 



^ tan ^ (y v/c/g) + = t . 



At t = 0 , y = 






Z ( 0 ) . There fore 

m 

X _i • 

D, = tan (Z (0 )/c7g) 
1 /7TTT m 



(D-22) 



(D-23) 



»/^ 



Solving for y and substituting for , 

y - /g/C tan(D^/CG) = -/g/C tan(t/CG) 



or 



y = -/g7c tan(t/^) + Z^(0) 
dZ 



(D-24) 



(D-25) 



Substituting y = and letting Z^(0) be constant, 



dZ = [-/gTc tan(t/^) + Z (0)]dt 

m 



and 



z = -/g/C 






log cos(t/^) + Z (0)t + , 



m 



(D-26) 



(D-27) 



Att=0, Z=Z (0). Therefore = Z (0) and 

m 2 m 

Z(t) = log^ cos(t/^) + Z^(0)t + Z^(0) 
C e mm 



(D-28) 



where Z (0) and Z (0) refer to the Z velocity and position 
mm 

of the missile at the time the prediction is made. 

Solving for the time when Z (t) is equal to the target 
elevation is done by substituting Z (t) = Z^ and using the 
Newton-Raphson method. 

The solution of the X and Y motion equations for impact 
point will be shown for the X equation only and the Y equa- 



tion will be written by inspection. 

mX = -K^X^ 



(D-29) 



dt 



cy‘ 



(D-30) 
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Ki 

where y = X and C = Separating variables and integrating 



TT- + D, = t . 
Cy 1 



(D-31) 



Att=0,y=X (0). Therefore 

m 



Dl = -- 



cx (0) 

m 



(D-32) 



dX 

Substituting y = ^ and letting be constant 



= Ct - CD 
dX ^1 * 



(D-33) 



Separating variables and integrating 



X = p log^(Ct - ) + D, 

X^(0) 

m 



(D-34) 



At t = 0 , X = 



X (0) . 
m 



There fore 



® X (0) 
m 



(D-35) 



and 

X(t) 



^ log^(Ct - i logg(- — ) + X^(0) (D-36) 



X^(0) 
m 



x^(0) 

m 



X(t) = i log^(Ct - T— i — )(-X„(0)) + X^(0) 
C X (0) “ ” 

m 



(D-37) 



or 



K tX (0) 

X(t) = ^ log [1 - -= — ] + X (0) 

K, ^e m m 



(D-38) 



For this method the prediction equations are 



K,tX (0) 

X = ^ log [1 — ] + X (0) 

p K, ^e m m 



(D-39) 



and 



^ K.tY„(0) 

"^ r , = ^ ” — 1 + . 

p •'Q mm 



(D-40) 
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APPENDIX E 



REFERENCE SYSTEMS FOR THE PROBLEM 



The reference systems for the problem are based on the 
position of the launch point of the aircraft involved. The 
inertial navigator of the aircraft is initialized when the 
aircraft is launched. All other points in the problem are 
also referred to this initialization point of the aircraft 
even though they may not be measured with respect to that 
point at first. The measured missile state vector is an 
example of this type in that it is first referred to the 
mother aircraft and then converted to a vector which is 
referred to the initialization point. 

The cartesian reference is a standard three-dimensional, 
right-hand reference system. Whenever a reference is shifted, 
the new frame is also a three-dimensional, right-hand system. 
Figure 22 illustrates the various reference points and systems. 

Azimuth angle is defined as the angle between true north, 
positive clockwise to the line-of-sight from the aircraft to 
the item of interest. In order to define the depression 
angle, it is first necessary to define the normal plane. The 
normal plane is defined as being a plane normal to a line pas- 
sing from the aircraft to the center of the earth. Depression 
angle is defined as the angle between the normal plane and the 
line-of-sight. This LOS must be contained in a plane which is 
perpendicular to the normal plane. These angles are shown in 
Fig. 23 and 24 respectively. 
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APPENDIX F 



THE MISSILE MODEL 



The missile model was relatively honest in its aerodynamic 
performance although it was not a true missile. The maximum 
yaw rate of the missile was limited to allow radial accelera- 
tion loads not to exceed ten G's. This was thought to be 
representative of current capabilities. 

The missile was powered by a rocket motor of 3000-lbs. 
thrust with a burn time of 30 seconds. After the burnout the 
missile glided for the remainder of its flight. 

The model was so designed that the mass and moment of 

inertia changed as the rocket motor burned. The total pro- 

/■ 

pellant weight was approximately 163.3 kg. The reduced mass 
and moment of inertia were then used for the remainder of the 
flight. A free body diagram of the missile appears in Fig. 

25. The equations of motion are 



mR = -D cosy - L siny + T cos 6 



(F-1) 



mZ = -D siny + L cosy + T sin6 - mG 



(F-2) 



I^^0 = -jL - jD sina + JICF - K^Q0 



(F-3) 



whe re 



Q = PV 



2 



(F-4) 



D = K^Q 



(F-5) 



L = K^Qa 



(F-6) 
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Figure 25. Missile Free-Body Diagram 



and 



T (thrust) 


13344. N 


L (lift) 


N 


D (drag) 


N 


V (velocity) 


m/sec 


CF (control force) 


N 


Y (velocity vector angle) 


rad 


a (angle of attack) 


rad 


0 (missile axis angle) 


rad 


• 

9 (axis rate) 


rad/sec 


9 (axis acceleration) 


rad/sec 


m (mass) 


453.6 kg 

(290.3 kg after burnout) 


G (gravity) 


9.8 m/sec^ 


i (control moment arm) 


2 . m 


j (center of pressure 

moment arm) 


1 . m 


I (moment of inertia) 

yy 


340.2 nt-m^ 

(217.7 nt-m^ after burnout) 


p (air density) 


1.284 kg/m^ 


(drag constant) 


.00556 nt/(kg/m^) (m/sec) ^ 


(control force 
constant) 

(damping constant) 


.01 nt/ (kg/m^) (m/sec) ^rad 

.004 nt/ (kg/m^) (m/sec) ^ (rad/sec) 


(lift constant) 


.008 nt/ (kg/m^) (m/sec) ^rad . 
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APPENDIX G 



THE LAUNCHING AIRCRAFT MODEL 

The launching aircraft model was a very simple model for 
the simulation. The aircraft was given an initial heading 
and velocity. It was then assumed to fly a straight line 
making no turns. The location of the aircraft was updated 
each sample period. 

Although preliminary investigation had shown that the 
most critical component of noise in the entire system was 
that in the inertial navigation package, the inertial navi- 
gator was considered to be noiseless. The reasons for 
neglecting the noise were that 1) the primary purpose of this 
investigation was in the field of application of impact pre- 
diction to an operator-controlled, air-to-surface weapon 
system, and 2) no figures were available concerning sample- 
to-sample noise in inertial navigators. 
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APPENDIX H 



MAN/SYSTEM INTERFACE 

In any system with a human operator, it is important that 
information be presented concisely and in such a manner that 
it can best be used by the operator. The communication link 
in any such application is a two-way channel. The system 
sends information to the operator in one direction and the 
operator-generated control commands travel in the other. 

The link to the operator carried two items of information. 
The first was the most recent prediction of impact point and 
the second was the predicted time until impact. 

The first item, impact point prediction, was displayed 
to the operator on a display scope with the target assumed to 
be in the center. This presentation was limited in such a 
manner that the predicted impact point could never appear 
off-scale. Figures 26, 27, and 28 give examples of the 
limiting action which was employed. The display, when in a 
limited state, gave no information as to the error of the 
predicted impact point except the direction in which it was 
off-scale . 

The display was generated such that the missile always 
approached the target from the bottom of the scope. This 
required a transformation of the axes of the display. This 
was necessary so that the motion of the control stick always 
drove the missile toward the target. 
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Figure 26. Impact Point Within Scale 
on Display 




Figure 27. Impact Point Off-Scale 
Left on Display 



f 


1 


1 





Figure 28. Impact Point Off-Scale 

Left, Off-Scale Vertically 
on Display 
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The transformation of axes consisted of both a translation 
and rotation. The translation was accomplished by generating 
a reference frame which had its origin located on the target. 
The predicted impact point was then referred to this origin 
and the axes rotated in such a manner as to make the Y axis 
of the rotated system coincide with the LOS from the current 
measured missile position to the target. This process is 
illustrated in Fig. 29 and 30. 

It can be shown that the location of the predicted im- 
pact point in the translated and rotated reference frame 
(display reference system) is 



X, = AX,^ 
— d — df 

where 



(H-1) 



and 




A 



cos 0 -sin 0 
sin 0 cos 0 




(H-2) 



(H-3) 



(H-4) 



0 = tan"^( (X. -X^)/(Y. -Y^) ) . (H-5) 

t m t m 

The second item of information presented to the operator 
was the predicted time remaining until impact. This value 
was computed once each sample period throughout the missile 
flight. It was generated on the digital computer, converted 
to an analog voltage, and displayed on a digital voltmeter. 
The display was scaled such that times were displayed only if 
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Figure 29. Translation of Axes for Display 




I 



Figure 30, Rotation of Axes for Display 
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they were less than or equal to 100 seconds. Any value of 
time which was greater than 100 seconds was displayed as 100 
seconds . 

Commands from the operator to the system formed the other 
direction of the two-way communication link. These commands 
were generated by a control stick which the operator used to 
guide the missile. Full deflections on the stick corresponded 
to full deflections on the missile controls. The pitch con- 
trol was scaled to provide 45 degrees of control movement on 
the missile when the stick was fully deflected. The yaw con- 
trol provided a yaw rate of 20 degrees/second when fully de- 
flected. This value was also subject to limiting by the 
missile dynamics. 
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APPENDIX I 



A DISCUSSION OF THE SIMULATION PROGRAM 

Tnis appendix describing the program used in the simula- 
tion of this problem is included in the event that some other 
investigator would desire to continue research in this oar- 
ticular field. In it the author has tried to include all 
pertinent information concerning the routines and their func- 
tions and especially the programming "tricks" which were used 
in some instances to provide a particular function in the 
overall program. 

The appendix is divided into four sections. The first 
section is a summary of the routines and the function which 
each performs. The second section lists the routines and 
the variable names within each routine. The method of inter- 
routine communication is the subject of the third section 
and the programming "tricks" used are discussed in a fourth 
section. A complete listing of the program and an analog 
patching diagram follow the appendix. 

A. SUMMARY OF THE ROUTINES AND THEIR FUNCTIONS 

1 . Main Line Program 

The main line program serves to perform two functions. 
The first is as a vehicle for input to the program and output 
from the program, and the second is as a calling and control- 
ling sequence for the various routines which do the actual 
computation for the simulation. If operated in a segmented 
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mode of core allocation, the main line program is always core 
resident . 

Several test functions are performed by the main line 
program also. It tests for the condition of problem termi- 
nation. Problem termination occurs when the measured Z of 
tne missile is equal to or less than the target Z. It also 
controls switching from the passive target location mode of 
the problem to the active mode in which the missile is fired 
and guided to the target. 

2 . Initialization Routines 

a. Inertial Navigator Initialization 

The subroutine INIT initializes the inertial 
navigation routine by setting the initial position of the 
aircraft including the altitude. It also initializes the 
heading and computes component velocities . The final items 
initialized by the routine are the elapsed time counter, 
missile flight time counter, and the problem time increment. 

b. Missile Initialization 

The subroutine MINIT initializes the missile model 
at the time of launch. All parameters of the missile model 
are set equal to the aircraft parameters at the time of launch. 
Also the routine gives initial values to the measured missile 
vector to enable other routines to function properly during 
the first second of the simulation until all routines are 
generating full outputs. 
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3. 



Aircraft Simulation Routine 



The subroutine NAVPCK is the routine which provides 
aircraft state vector information for the' problem. The 
routine provides a straight-line flight path for the aircraft 
at any initial heading from 0 to 359 degrees at any initial 
velocity. The routine receives flight path information in 
R/0 form and provides flight-path information in the three 
cartesian coordinates. 

4 . Passive Target Location 

The subroutine PFIL simulates the action of the hybrid 
filter presented in section II-A-3. This location action is 
simulated by interpolating and normalizing values of variance 
provided for the hybrid method and generating noise which is 
added to the true target location. The final noisy location 
of the target is then used throughout the remainder of the 
problem. 

5 . Missile Model 

The missile model routine is of sufficient importance 
that an entire previous appendix was devoted to it. The reader 
is referred to Appendix F for a discussion of the missile model. 

6 . Routines Involved in Missile Tracking 
a. Radar Vector Generation 

The subroutine REVCON is the routine which gener- 
ates the true radar vector from the geometrical relationships 
between the true missile and the aircraft. This vector rep- 
resents the true radar state before any noise is added and 
before quantization takes place. 
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b. Noise Addition and Quantization of Radar Vector 
The subroutine RADAR2 takes the true radar vector 

produced by REVCON and performs two operations on it. First 
the routine adds the measurement noise associated with the 
particular measured state variable. Second, the routine 
quantizes the measurements into the discrete form which 
would be provided by a digitized radar output unit. 

c. Radar Vector Filtering 

The subroutine AFIL provides the filtering of 
tne noisy, quantized radar vector to produce a radar vector 
whose variance is less than that of the noisy vector. The 
subroutine uses a schedule of precomputed filter gains and 
a switch which determines the number of samples after which 
the gain schedule is trUhCated. This value can be from one 
to thirty seconds assuming gains computed for a sample time 
of one second. 

d. Spherical tO Cartesian Coordinate Conversion 
The subroutine CORDCN processes the filtered 

radar vector spherical values to produce a measured missile 
vector composed of missile position and velocities relative 
to the true problem reference. The subroutine first pro- 
duces a cartesian vector relative to the aircraft. The air- 
craft component positions and velocities are then added to 
the relative components to produce the final measured missile 
vector. 
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7. 



Initial Missile Control Routine 



The subroutine SHOOT controls the flight of the mis- 
sile during the boosted phase (first thirty seconds) of its 
flight. The routine attempts to drive the missile velocity 
vector to assume an up angle of 47.5 degrees and to turn the 
missile velocity vector projection on the XY plane to coin- 
cide with the LOS from the missile to the target. 

8 . Impact Prediction and Display 

a. Impact Point Prediction 

The subroutine PREDCT develops the X and Y im- 
pact point prediction values. This routine employs the sub- 
routine NR to assist in this task. 

b. Time-Until-Impact Solution 

The subroutine NR solves for the time to go 
until impact for the missile by using the Newton-Raphson 
iterative method of solution for differential equations, 

c. Command Generation 

The subroutine COMGEN converts the outputs of 
the control stick to values which can drive the missile mod- 
el. This routine uses the SDS 9300 routines for analog-to- 
digital conversion. Full-scale values of pitch correspond 
to plus or minus 45 degrees of control deflection and full- 
scale values of yaw correspond to turn rates of plus or 
minus 20 degrees/second. 

9 . Prediction Display and Statistics Routines 

a. Display Generation 

The subroutine DISGEN converts the predicted im- 
pact points into values which can be used to drive the 
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situation display. This routine performs the axis rotation 
and translation necessary to generate a useful display. The 
limiting of the display range is also performed, 
b. Minimal Statistics Routine 

The subroutine MISS provides minimal statistics 
about the success of the run at the time of flight termina- 
tion. This routine provides the X and Y miss distances, 
impact time of the missile, and total miss distance. More 
complete statistics are generated by a separate program which 
analyzes data which is written on a storage medium during the 
execution of the program. The separate analysis routine pro- 
vides such results as missile flight profile, control applied 
vs. time, missile location errors vs. time, and others. 



B. ROUTINE SYMBOLS AND THEIR MEANINGS 



1 . Main Line Program 
Symbol 
GBQ/KBQ 

RUNS/RNS 

NP 

SMU, ARNG, ARNGD, AAZ , 
VL, THI, Z, TME 

ETIME, TIME, FTIME 
ETA 



Meaning 

Time for active filter 
gain truncation. 

Number of runs for the 
program . 

Switch used to initialize 
AFIL. 

ADP Mean of noise and noise 
variances . 

Initial velocity, heading 
(in degrees) , aircraft 
altitude, and time incre- 
ment. 

Elapsed time, time incre- 
ment, flight time. 

Initial launch angle. 
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2 . 



Subroutine INIT 



Symbol 


Meaning 


XA , YA , Z A , XDA , YDA , Z DA 


Aircraft location and 
velocity in cartesian 
coordinate frame. 


CK,CT,TH,V 

3. Subroutine NAVPCK 


Conversion factor (degrees 
to radians) , converted 
heading angle, heading in 
degrees , velocity in 
meters/ second. 


Symbol 


Meaning 


CK ,CT 


Conversion factor (degrees 
to radians) , converted 
heading . 


XA,YA,ZA , XDA , YDA , ZDA , TH , V 
4. Subroutine MINIT 


Aircraft state vector. 


Symbol 


Meaning 


YAW; PITCH 


Missile command vector. 


XMM ; YMM , ZMM ; XDMM ; YDMM ; 
ZDMM 


Measured missile vector 
in position and velocity. 


XM ; YM , ZM , XDM , YDM ; ZDM ; 
XDDM,YDDM,ZDDM; 
HEAD ; GAMMA , THET ; 
THETD;THETDD 

5. Subroutine PFIL 


True missile vector. 


Symbol 


Meaning 


X,Y 


True target location. 


QK,QM,K,M,J 


Constants used for 
interpolation . 


XS ;YS ;A 


X variance, Y variance, 
noise value returned by 
GAUSS . 


XT,YT;CA 


Noisy target location, 
CPA in problem. 
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Subroutine MISAD 



6 . 



Symbol 


Meaning 


RHO,T ,G 


Air density, thrust, 
gravity . 


RE FAR 


Aerodynamic reference 
area . 


K1 ,K3 ,K4 ,K5 


Drag, lift and moment 
constants . 


K , YINC , ANG , VEL , VELSQ 


Advance count, yaw incre- 
ment, current missile 
heading, current velocity, 
current velocity squared. 


DELTA, GAMMA, ALPHA 


Control angle, missile 
axis angle, angle of 
attack . 


D,QL,QM 


Drag, lift, and moment 
forces . 


QMASS,QIYY,QMULT 
7. Subroutine REVCON 


Missile mass, moment of 
inertia, constant relating 
mass and moment. 


Symbol 


Meaning 


X,Y,Z 


Work variables. 


RPRI ,RDPRI ,T1 ,T2 ,VLXY 


Work variables for develop 
ing the radar vector. 


R,RD,T ,P 

8. Subroutine RADAR2 


True radar vector. 


Symbol 


Meaning 


Cl 


Bin size for 13-bit shaft 
encoder . 


I TEMP ,AN 


Temporary variable used 
in the conversions , and 
noise value returned by 
GAUSS . 


R , RD , T , P 


Noisy, quantized radar 
vector . 
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Subroutine AFIL 



9 . 



Symbol 


Meaning 


NP 


Switch to control filter 
initialization. 


ZRl ,ZR2 ,ZA1 ,ZD1 


Measurements of range, 
range rate, azimuth, and 
depression angles. 


XRl , XR2 , XAl , XA2 , XDl , XD2 


Estimated values of range, 
range rate, azimuth, 
azimuth rate, depression, 
and depression rate. 


GRG^GAZ ,GDP 

10. Subroutine CORDON 


Range , azimuth , and de- 
pression filter gains. 


Symbol 


Meaning 


CSP,SNP,CST,SNT 


Sine and cosine of (j) 
and 0 . 


CSPSNT , CSPCST ,SNPSNT , 
SNPCST 


Cross products of sine 
and cosine values. 


R,RD,T,TD,P ,PD 


Range, range rate, theta, 
theta rate, phi, and phi 
rate . 


XMM , YMM , ZMM 


Measured missile cartesian 
position . 


XDMM,YDMM,ZDMM 
11. Subroutine SHOOT 


Measured missile cartesian 
velocity . 


Symbol 


Meaning 


GAMMA, HEAD 


Current missile velocity 
vector elevation and cur- 
rent missile heading. 


ETA,A1 


Desired missile velocity 
vector elevation and cur- 
rent LOS angle from the 
missile to the target. 


YAW, PITCH 


Missile command vector. 
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12. Subroutine PREDCT 





Symbol 


Meaning 




QMASS ,RHO,REFAR 


Missile mass, air density/ 
missile aerodynamic refer- 
ence area. 




K2 ,K3, V 


Drag coefficients and 
missile velocity. 




XP ,YP , ZP 


X and Y impact position 
predictions, and time 
until impact. 


3. 


Subroutine NR 






Symbol 


Meaning 




QMAS S , RHO , REFAR , G , V 


Missile mass, air density, 
missile aerodynamic refer- 
ence area, gravity, and 
missile velocity. 




K1 , ALPHA 


Drag coefficient and 
exponential argument. 




A,B ,C,TOLD,TNEW 


Coefficients of prediction 
equations , last time esti- 
mate , new time estimate. 




FT , FDT 


Function value , and 
derivative . 


14. 


Subroutine COMGEN 






Symbol 


Me aning 




YAW, PITCH 


Missile Command vector. 




Numeric arguments 


Conversion factors. 


15. 


Subroutine DISGEN 






Symbol 


Meaning 




XDF,YDF 


Translated predicted 
impact point. 




XTA,YTA 


Used to determine 
rotation angle. 




ALPH 


Rotation angle. 
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XD,YD,ZD 



16 . Subroutine MISS 
Symbol 

ZTMl ,TI1 ,TI2 
XIMl ,YIM1 ,XI ,YI 



ETI,FTI 
K1 ,K2 ,QK2 

XMS ,YMS ,XT ,YT 



Display points and time 
until impact. 



Meaning 

Last preimpact missile Z 
position, times which 
bracket impact time. 

Last preimpact X and Y 
positions, and missile 
impact positions in X 
and Y . 

Impact time as elapsed 
time and flight time. 

Variables used to edit 
inaccuracies out of 
flight time clock. 

X and Y miss distances, 
and last estimated target 
position . 



C. INTERROUTINE COMMUNICATION 

Interroutine communication is carried out throughout th« 
program by means of labeled FORTRAN COMMON blocks. These 
blocks are designed in such a manner that variables relatina 
closely to each other have been grouped into sections named 
to reflect their functions. The block names are given below 
along with the variables included and a description of the 
function of the block. 

1. Blank COMiMON 

The variables included in blank COMMON are IX, ETIME, 
TIME, and FTIME. IX was originally the seed value for the 
GAUSS routine which produces random numbers. It was not re- 
tained when a change in random number generators was made. 
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ETIME , TIME, and FTIME are the problem elapsed time, time 
increment, and missile flight time respectively, 

2 . Block TGPRD 

The variables included in the block TGPRD are the 
cartesian location of the target and the X and Y predicted 
impact points of the missile. ZP is the time until impact. 

3. Block SIGMAS 

This block contains the values of the sigmas to be 
used in the radar to introduce noise in the measurements. 
These values are set in the main line program. 

4 . Block GAINS 

The computed filter gains used by the active radar 
filter are stored in this dimensioned block. The block 
provides storage for 120 range filter gains and 60 azimuth 
and depression filter gains each. 

5. Block MISLM 

This block contains the values of measured missile 
position and velocity developed by the missile tracking 
routine . 

6 . Block ACRFT 

The values of aircraft position and velocity at each 
time k are held in this block. Also this block contains the 
aircraft velocity in meters per second and aircraft heading 
in degrees . 

7. Block MISLT 

This vector contains all values which pertain to the 
true missile. The vector contains position, velocity, and 
acceleration in cartesian coordinates; missile axis angle. 
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angular velocity, and angular acceleration; and velocity vec- 
tor heading and elevation angle. This is updated each sample 
period to reflect the true state of the missile. 

8. Block COMND 

The commands to the missile in yaw and pitch comprise 
the two elements of this block. The yaw values are up to 
twenty degrees per second in either direction and pitch val- 
ues range from plus to minus forty -five degrees in control 
deflection . 

9 . Block SIGS 

The values of variance which are used in the passive 
location routine PFIL are transferred to the routine by this 
block. Each of the values is dimensioned to accept fifty- 
nine values of variance. 

10. Block EADAR 

This block is used to transfer the values of the radar 
information through all of the routines of the missile track- 
ing portion of the program. The six elements of the vector 
are range, range rate, azimuth, azimuth rate, depression, and 
depression rate. 

D. PROGRAMMING "TRICKS" 

Only one programming irregularity exists in the simula- 
tion program. In the subroutine PFIL the variable ZP in the 
COMMON block TGPRD is used as a switch to allow the routine 
to hold the true target location after the first call. The 
variable ZP was chosen since it was already in the COMMON 
block and since it is not used by any other routine in the 
passive phase of the program. 
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PAUSE 77 

61 IF(TEST(7)) 60/6li6l 
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